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ABSTRACT 


As  a  back-up  navigation  method  for  Polaris  submarines,  the  feasi¬ 
bility  of  moon  doppler  navigation  on  a  bistatic  basis  is  investigated.  The 
problem  is  that  of  finding  the  latitude  and  longitude  of  a  vessel  by  meas¬ 
uring  doppler  shift  and  time  rate  of  change  of  doppler  shift  in  radio 
waves  transmitted  from  a  fixed  site  and  reflected  from  the  moon. 

Fora  prescribed  navigation  accuracy,  one  must  determine  the 
necessary  precision  in  measuring  doppler  frequency  and  time  rate  of 
change  of  doppler  frequency.  A  sample  calculation  is  made  for  a  location 
on  the  western  edge  of  the  Norwegian  Sea,  and  it  is  found  that  if  the 
probable  error  in  latitude  is  specified  to  be  0.8  naut  mi  and  the  probable 
error  in  longitude  to  be  1.3  naut  mi,  the  probable  error  which  can  be 
allowed  in  doppler  shift  is  0.796  cps  and  the  permissible  probable  error 
in  time  rate  of  change  of  doppler  shift  is  9.511  x  10”5  cpsps,  with  a 
transmitted  frequency  of  2200  Me. 

The  calculated  value  for  doppler  frequency  was  2459  cps,  and  that 
for  doppler  rate  was  0.1381  cpsps.  This  means  that  doppler  frequency 
must  be  measured  to  an  accuracy  of  0.0324  percent  and  doppler  rate  to 
an  accuracy  of  0.0689  percent,  for  this  particular  example.  The  meas¬ 
urement  of  doppler  frequency  seems  within  the  realm  of  technical 
capability  in  the  near  future.  The  accuracy  of  the  measurement  of  dop¬ 
pler  rate  is  still  in  question.  However,  it  appears  that  the  rate  can  be 
computed  within  the  required  accuracy  by  fitting  an  approximate  ana¬ 
lytical  expression  to  the  set  of  points  representing  the  measured  doppler 
shift  and  differentiating  the  expression,  or  by  using  finite  difference 
methods  with  the  observed  data  to  get  the  time  derivative. 

Propagation  effects  due  to  variation  of  index  of  refraction  in  the 
atmosphere  have  not  been  taken  into  account.  These  effects  may  be 
important,  and  should  probably  be  the  subject  of  investigation  in  the 
near  future. 
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continuing. 
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NAVIGATION  BY  MOON  DOPPLER  EFFECT 


INTRODUCTION 

It  is  desirable  that  a  back-up  navigation  system  be  provided  for  the  Polaris  submarine, 
to  supplement  the  basic  ship  inertial  navigation  system  (SENS).  In  this  report  equations  are 
derived  which  show  how  the  latitude  and  longitude  of  a  vessel  can  be  obtained  from  measure¬ 
ments  of  doppler  shift  and  time  rate  of  change  of  doppler  shift  in  electromagnetic  waves  ema¬ 
nating  from  a  source  at  a  known  fixed  location  and  reflected  from  the  moon.  Further  equations 
are  derived  which  show  the  accuracy  of  measurement  of  these  two  quantities  required  to 
obtain  prescribed  accuracies  in  latitude  and  longitude. 


DOPPLER  EQUATIONS 

The  classical  equations  for  the  doppler  shift  for  a  moving  source  and  for  a  moving 
observer  are  as  follows  (1): 

1.  For  a  moving  source 


f 


R 


(l) 


where 


fT  =  radiated  frequency 

fR  =  frequency  received  by  observer 

v T  =  velocity  of  source 

c  =  velocity  of  propagation  of  radiation.* 

Here  vT  is  taken  as  positive  when  the  source  moves  toward  the  observer,  since  the  observer 
then  receives  a  higher  frequency. 

2.  For  a  moving  observer 


where  vR  =  velocity  of  observer.  Here  vR  is  taken  as  positive  when  the  observer  is  moving 
toward  the  source. 


Changes  in  velocity  of  propagation  through  the  atmosphere  have  not  been  taken  into  account. 


1 
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3.  Consequently,  when  both  Observer  and  source  are  in  motion, 


Since  c  is  very  large  with  respect  to  the  other  velocities  involved, 


fR  =  fT 


c 


(3) 


(4) 


Suppose  that  upon  the  earth,  considered  stationary,  a  transmitter  radiates  a  signal 
of  frequency  fT  and  this  signal  is  reflected  by  the  moon.  If  the  moon  is  moving  at  a  speed 
vT  with  respect  to  the  transmitter  site,  the  reflected  signal  will  have  a  frequency 


F 


M 


(5) 


If  the  reflected  signal  is  received  at  some  point  on  the  earth,  the  moon  may  now  be  con¬ 


sidered  to  be  a  source,  in  motion  toward  the 
have 


observer,  and  radiating  at  frequency  fm.  We 


where  FR  is  the  received  frequency  and  vR  is  the  speed  of  the  moon  with  respect  to  the 
receiver.  Then  the  doppler  shift  is 

Af  =  FR  -  fT.  (8) 

Equations  (7)  and  (8)  combined  give 


~  fT 

Af  =  —  (VR  +  VT)  ' 


(9) 


MOTION  OF  A  POINT  ON  EARTH  RELATIVE  TO  THE  MOON 

Trotter  (2)  has  derived  a  good  approximation  for  the  motion  of  a  point  on  the  earth 
relative  to  the  moon.  This  equation  is,  for  the  receiving  site, 
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where 


vR  =  co  R  cos  8  cos  Lr  sin  HR 

co  =  angular  velocity  of  the  earth's  rotation 
R  =  equatorial  radius  of  the  earth 
8  =  declination  of  the  moon 
Lr  =  latitude  of  receiving  site 
Hr  =  hour  angle  of  receiving  site. 


Similarly,  for  the  transmitting  site, 

vT  =  co  R  cos  8  cos  Lt  sin  Rj. 


(10) 


(11) 


where 


Lt  =  latitude  of  transmitting  site 
Ht  =  hour  angle  of  transmitting  site. 

The  purpose  of  employing  these  expressions  for  vR  and  vT  is  to  obtain  Af  by  substitut¬ 
ing  in  Eq.  (9).  The  transmitted  frequency  fTis  known.  The  value  of  co  is  computed  from 
the  length  of  the  sidereal  day,  which  is  given  by  Smart  (3)  as  23h  56m  4s. 091. 

Smart  may  also  be  consulted  for  the  definition  of  the  astronomical  terms  used.  The 
earth's  radius  and  the  coordinates  of  the  transmitting  site  are  taken  from  Ref.  2  and  repro¬ 
duced  in  Appendix  A.  In  addition,  this  reference  gives  values  and  probable  errors  for  all 
the  quantities  involved,  except  LR  and  hr.  The  declination  8  for  any  date  and  time  can  be 
obtained  by  consulting  the  American  Nautical  Almanac. 

Strictly  speaking,  the  value  of  co  used  should  be  the  difference  between  for  the  earth 
and  &  for  the  moon,  which  we  can  call  <^>E  and  a ^ .  The  value  of  ojm  for  a  circular  lunar 
orbit  is  2.661697  x  10~e  rad/sec,  which  is  only  3.56  percent  of  .  The  variation  in  ojm  due 
to  the  ellipticity  of  the  lunar  orbit  amounts  to  10.9  percent  of  Hence  varies  from 
3.25  to  4.05  percent  of  oj£.  Therefore,  is  neglected  (see  Appendix  B). 


CALCULATION  OF  LATITUDE  AND  LONGITUDE  OF  A  VESSEL 

Let  us  proceed  to  derive  expressions  for  LR  and  hr  based  on  measurements  of  doppler 
shift  and  time  rate  of  change  of  doppler  shift.  For  convenience,  we  group  the  necessary 
previous  equations 


vt) 


(9) 


vR  =  co  R  cos  8  cos  Lj^  sin  HR 


(10) 


V 


T 


co  R  cos  8  cos  Lj.  sin  Rp. 


(ii) 
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Then 


and 


fT  a>  R  cos  S 


Af  = 

c 

-(cos  Lp  sin  Hp  +  cos  LT  sin  H^J 

(12) 

*T  R  cos 

8  . 

(13) 

Af  = 

c 

“  (cos  Lp  cos  Hp  +  cos  Lt  cos  Rj. 

The  extra  w  in  Eq.  (13)  appears  because 


dHp  dKj. 
dt  “  dt 


Let 


A  = 


c  A  f 

fT  CO  R  COS  8 


B  -  cos  Lt  sin  Ep 
c  Af 

D  =  - 

fT  a >2  R  cos  8 

E  =  cos  Lt  cos  Rj. . 

Then  Eqs.  (12)  and  (13)  can  be  written 


A  -  B  -  cos  Lr  sin  HR 


D  -  E  =  cos  LR  cos  Hp . 


(14) 

(15) 

(16) 

(17) 

(18) ' 

(19) 

(20) 


Hence 

A  -  B 

tan  =  r~T 

D  -  E 


co  c  A  f  -  fx  co2  R  cos  8  cos  L_  sin  H~ 

= - : - - - - - -X  (21) 

c  Af  -  fT  co2,  R  cos  8  cos  LT  cos  Hr 

The  expression  for  cos  Lp  could  be  obtained  by  squaring  Eqs.  (19)  and  (20)  and  adding, 
but  a  simpler  expression  results  when  Eq.  (2)  is  used  to  get  LR  in  terms  of  HR : 


D  -  E 

cos  Hr 


cos  1^  = 


(22) 
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When  cos  hr  and  D  -  e  are  negative  (as  in  the  numerical  example  used  later)  Eq.  (22) 
can  be  written 


.  _  E  "  D 

cos  Lp  -  - . 

|  cos  Hjjl 


Using  the  identity 


cos  Lr  =  ^1  +  tan2  1/2 


(23) 


(24) 


and  Eqs.  (17)  and  (18),  Eq.  (23)  becomes 

cos  Lp  =  ^1  +  tan2 


hr)  1/2  (c 


cos  Lt  cos  Rj. 


C  Af 


fT  co 2  R  cos  8 


;) 


(25) 


The  relation  between  the  longitude  and  the  local  hour  angle  HR  of  the  moon  is 

GHA  -  Hr  =  longitude  (26) 


where  GHA  is  the  Greenwich  hour  angle  of  the  moon,  obtainable  from  the  American  Nautical 
Almanac  for  any  particular  date  and  time.  Therefore,  if  we  know  hr,  we  know  the  longitude 
of  the  vessel. 

Since  the  tangent  and  the  cosine  are  periodic  functions,  the  results  given  by  Eqs.  (21) 
and  (25)  are  not  unique  in  the  purely  mathematical  sense.  However,  since  it  may  be 
expected  that  latitude  is  known  within  90  degrees  and  longitude  within  180  degrees,  the 
answers  are  unique  for  navigational  purposes.  Equations  (21)  and  (25)  are  quite  simple 
in  form;  hence,  it  should  not  be  difficult  to  design  a  computer  to  obtain  specific  numerical 
results. 


PROBABLE  ERROR  ANALYSIS 

It  is  now  necessary  to  make  a  probable  error  analysis  to  find  out  how  accurate  doppler 
and  doppler-rate  measurements  have  to  be  to  provide  a  specified  accuracy  in  latitude  and 
longitude.  From  Eq.  (12), 

fT  co  R  cos  8  /  \ 

Af  =  - - -  (cos  sin  HR  +  cos  LT  sin  HT j  (27) 


and,  from  Eq.  (13), 


f~,  co2  R  cos  8 


cos  Lr  cos  Hr  +  cos 


Hr)  * 


Af  = 


C 


(28) 
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Appendix  A  lists  the  values  of  the  parameters  and  their  probable  errors. 


NUMERICAL  EXAMPLE 

For  the  purpose  of  carrying  out  a  sample  computation,  let  us  arbitrarily  choose  the 
location  of  the  receiving  site  (Polaris  vessel)  at  the  western  edge  of  the  Norwegian  Sea, 
at  latitude  LR  =  70°N  and  longitude  =  10°  W.  The  longitude  is  considered  positive  if  the 
point  considered  is  west  of  the  Greenwich  meridian.  Thus  we  may  write  the  longitude 
mentioned  above  as  +10°.  The  relation  between  the  longitude  and  the  local  hour  angle 
H„  of  the  moon  is,  as  mentioned  before, 

GHA  "  HR  =  longitude.  (26) 

For  the  sample  calculation  Trotter’s  data  for  July  7,  1956,  at  1900  hours  Universal 
Time  are  used.  At  this  time  GHA  -  110°  11'.  1  (see  Ref.  2b  or  consult  the  American  Nautical 
Almanac).  Hence  we  have  the  equation 

no°ii'  .1  -  hr  =  +10°  (31) 

or 

Hr  =  100°ll'.l.  (32) 


=  100°ll'.l. 


(33) 
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Let  us  also  assume  a  transmitted  frequency  fT  =  2200  Me.  Basically,  it  is  necessary 
to  get  the  permissible  probable  error*  in  Af  and  A-f  corresponding  to  a  prescribed  probable 
error  in  latitude  and  longitude  (or  hour  angle).  It  is  also  of  interest  to  find  the  actual  values 
of  Af  and  Af  .  Using  the  parameter  values  given  in  Appendix  A  and  substituting  in  Eqs.  (27) 
and  (28),  we  get 

Af  =  2459  cps  (34) 


and 


Af  =  0. 1381  cpsps . 


(35) 


For  the  error  analysis,  we  use  Eqs.  (29)  and  (30)  and  get  (see  Appendix  C) 


and 


fiAf  =  0.796  cps 


(36) 


SA-f  =  9.511  x  10"5  cpsps  (37) 

assuming  t,  =  10"10f_,  ?,  =  0.8  nautical  miles,  and  S  =  1.3  nautical  miles.t 

*  T  T  R  HR 

Dividing  8Af  by  Af ,  we  get  3.24  x  10“4;  dividing  8^f  by  Af,  we  get  6.89  x  10“4.  This 
means  thatdoppler  frequency  must  be  measured  to  an  accuracy  of  0.0324  percent  and 
doppler  rate  to  0.0689  percent. 


The  measurement  of  doppler  frequency  seems  to  be  within  the  realm  of  technical 
capability  in  the  reasonably  near  future.  The  accuracy  of  doppler- rate  measurement  is 
still  in  question.  However,  it  appears  that  the  rate  can  be  computed  within  the  required 
accuracy  by  fitting  an  approximate  analytical  expression  to  the  set  of  points  representing 
the  measured  doppler  shift  data  and  differentiating  the  expression,  or  by  using  finite  dif¬ 
ference  methods  with  the  observed  data  to  get  the  desired  time  derivative.  This  conclusion 
follows  from  the  following  analysis:  If  we  have  two  quantities  x  and  y  with  probable  errors 
8  and  g  ,  the  maximum  probable  error  in  the  quotient  x/y  will  be 

x  y 


y 

X 


(38) 


Let  us  select  two  values  of  time,  tj  and  t2  ,  with  corresponding  doppler  frequencies 
Af  x  and  Af  2 .  We  can  then  form  the  difference  quotient 


Af j  -  Afa 
~1  ~ ^ 


Many  people  prefer  to  use  the  standard  deviation  cr .  The  relation  between  the  two  is 


Probable  Error  =  0.6745  cr. 

^Strictly  speaking,  latitude  and  longitude  are  angles.  When  they  are  expressed  in  terms 
of  linear  measure  the  distances  given  are  those  which  subtend  the  actual  angular 
latitude  and  longitude  at  the  center  of  the  earth. 
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and,  from  Eq.  (38), 


/Afj  -  Af\  _  Afi  -  Af 2 

8(Af,  -  Af2)  £( t J  -  t2) 

....  .  .  + 

V*1  ”  t2 j  tl  ~  1 2 

_  ^  1  ^2  t  1  ”  t2  _ 

(39) 


(Here  the  argument  notation  is  used  instead  of  the  subscript  notation  for  convenience  in 
printing.) 

If  the  permissible  probable  error  in  Af  is,  as  given  above,  3.24  x  10"4  x  Af ,  then,  since 
we  are  considering  two  doppler  frequencies  at  two  different  times,  the  maximum  probable 
error  in  their  difference  will  be  2  x  3.24  x  10‘4  (Afx  -  Af2)  =  6.48  x  10'4(Af1  -  Af2)  . 

As  for  errors  in  time  measurement,  let  us  assume  (4)  that  we  will  be  able  to  measure 
a  given  frequency  to  an  accuracy  of  1  part  in  10l°.  That  is  to  say,  let  us  take 

S{  =  10_10f  (40) 

where  f  -  frequency  in  cps  and£f  =  probable  error  in  frequency  in  cps.  At  frequency  f  cps, 
one  cycle  is  f-1  =  t  seconds  long.  Then,  by  the  previously  used  error  analysis  method, 


or 

et  =  r2  s{  =  f-2  x  io~10  f  =  io^V1  =  io"10t.  (42) 


This  is  to  say  that  the  probable  error  in  frequency,  divided  by  frequency,  is  the  same  as 
the  probable  error  in  time,  divided  by  time.  Hence  the  maximum  value  of  the  probable 
error  S(tA  -  t2)  in  Eq.  (39)  is  2  x  10"lo(t1  -  t2)  .  If  the  time  interval  is  sufficiently 
small,  then  the  difference  quotient  is  approximately  equal  to  the  derivative: 

p.  =  Af  (6.48  X  10"4  +  2  X  1(T10).  (43^ 

cAf 

The  second  term  in  the  parentheses  may  be  considered  negligible  in  comparison  with  the 
first.  Then,  since  Af  =  0.1381  cpsps  the  calculated  error  is 

SAf  =  8.95  x  10”5  cpsps.  (44) 

Since  the  prescribed  error  is  of  the  same  order  (SA*f  =  9.511  x  10"5  cpsps)  it  would  seem 
that  the  desired  accuracies  are  attainable.  Curves  have  been  plotted  showing  how  8Af  and 
SA*f  vary  as  SLr  and  8Hr  are  varied  (see  Appendix  D,  Figs.  D1  -  D4). 

It  is  of  interest  to  calculate  the  residual  errors  in  latitude  and  hour  angle  when  we 
assume  that  SAf  and  SA*f  are  zero.  Using  the  set  of  parameters  previously  assumed,  we 
obtain  (see  Appendix  E) 


^  ^  -  0.616  naut  mi 

(45) 

(Ew  ”  0.277  naut  mi. 

\  R/o 

(46) 

where  the  zero  subscripts  indicate  that  these  are  residual  errors. 
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These  results  show  that  there  would  be  a  finite  probable  error  in  latitude  and  longitude, 
even  if  Af  and  Af  were  to  be  measured  exactly;  this  is  because  of  unavoidable  errors  in 
measuring  other  physical  quantities  such  as  c,  R,  and  so  forth.  The  calculated  errors 
shown  in  Eqs.  (45)  and  (46),  however,  are  smaller  than  the  navigation  error  requirement, 
and  therefore  will  not  be  limiting  errors. 

CONCLUSIONS 

It  appears,  on  the  basis  of  the  analysis  involved  in  this  report,  that  moon  doppler 
navigation  is  feasible.  The  required  accuracy  of  doppler  frequency  measurement  is 
probably  attainable  in  the  near  future.  By  curve-fitting  or  finite  difference  methods,  it 
is  believed  that  doppler  rate  can  be  calculated  within  the  required  accuracy.  It  must  be 
noted  that  there  has  been  no  consideration  in  this  report  of  propagation  variation  in  the 
atmosphere. 

There  are  two  approaches  which  could  be  considered  in  an  attempt  to  make  the  method 
more  immediately  feasible.  One  is  to  relax  the  requirements  on  navigational  accuracy. 
According  to  the  curves  in  Appendix  D,  a  large  relaxation  in  required  navigation  accuracy 
would  result  in  relatively  small  increases  in  allowable  doppler  and  doppler  rate  error. 

The  other  way  of  improving  the  method  is  to  improve  electronic  techniques  to  the  point 
where  the  required  accuracies  can  be  more  reliably  attained.  It  appears  that  this  is  the 
direction  in  which  further  efforts  should  be  made. 
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Appendix  A 

LIST  OF  PARAMETERS  AND  THEIR  PROBABLE  ERRORS* 

c  =  2.997893  x  108 meters/second  (velocity  of  propagation) 

£c  =  5  x  102  meters/second  (probable  error  in  c) 

co  =  7.292116  x  10'5 rad/sec  (angular  velocity  of  earth  with  respect  to  moon) 

£  =  1  x  1CT8 rad/sec  (probable  error  in 
R  =  6.378260  x  106  meters  (equatorial  radius  of  earth) 

£  =  1  x  102 meters  (probable  error  in  r) 

§  =  20°  19 '(declination  of  moon,  July  7,  1956,  at  1900  UT) 

£s  =  4  x  10"7  radians  (probable  error  ins) 
sin  5  =  0.34721 
cos  8  =  0.93779 

Lt  =  38.545  degrees  =  0.673  radians  (latitude  of  transmitter) 

=  5  x  10”7  radians  (probable  error  inLT) 
sin  Lt  =  0.62303 
cos  Lt  =  0.78212 

Hj.  =  31.914  degrees  =  0.557  radians  (local  hour  angle  of  moon  at  transmitter) 

F  =  1  x  10-6  radians  (probable  error  inHT) 
ht 

sin  Rj.  =  0.55194 
cos  Rj,  =  0.83389 
Lr  =  70  degrees 

P  =0.8  naut  mi  =  2.323  x  10'4  radians^ 

lr 

sin  Lr  =  0.93969 

*Most  of  these  values  are  taken  from  B.  E.  Trotter,  “Navigation  Systems,”  NRL  Report 
5005  (Secret  Report,  Unclassified  Title),  pp.  19  and  20,  Oct.  22,  1957 

^The  conversion  of  nautical  miles  to  radians  is  derived  in  Appendix  C. 
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cos  Lp  =  0.34202 
Hp  =  100°  ll7. 1 

eu  =1.3  naut  mi  =  1.104  x  10”  3  radians^ 
hr 

sin  Hp  =  0.98425 
cos  HR  —  -0.17683 
fT  =  2200  Me 

=  2.200  x  10#  cps 
Sf  =  10-*°fT  =  0.22  cps. 


t 


The  conversion  of  nautical  miles  to  radians  is  derived  in  Appendix 


C. 


Appendix  B 

VARIATION  IN  ANGULAR  VELOCITY  OF  MOON 


The  angular  velocity  ^  of  the  moon  is  proportional  to  1/r2,  where  r  is  the  radius  of 
the  orbit  at  a  particular  time.  As  given  by  Smart,*  the  eccentricity  e  of  the  moon's  orbit 
is  0.0549.  From  the  same  sourcet  we  find  that  the  mean  geocentric  radius  of  the  orbit, 
which  is  the  same  as  the  semimajor  axis  a,  is  3.84400  x  108km. 

Using  the  relation 


a(l  +  e)  >  r  >  a(l  -  e)  ,  (Bl) 

and  the  given  figures,  it  is  found  that  the  mean  value  of  is  2.661697  x  10"*  rad/sec, 
which  is  3.6501  percent  of  the  earth’s  angular  velocity  o;E.  The  value  of  varies  by  10.9 
percent  from  the  mean  value,  which  means  that  its  value  is  3.6501  ±0.3979  percent  of  the 
earth’s  angular  velocity.  That  is  to  say,  varies  from  3.25  to  4.05  percent  of  oje  . 


W.  M.  Smart,  “Foundations  of  Astronomy,”  Longmans,  Green  and  Co.,  p.  176,  1947 
(4th  impression) 

V  258 
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Appendix  C 


DERIVATION  OF  8Af  AND  8Af 


The  doppler  shift  is 

f  _,Cc>R 

Af  =  “  cos  8  (cos  1^  sin  HR  +  cos  LT  sin  Rj.)  (Cl) 

For  the  parameters  assumed  (Appendix  A),  we  get 

Af  =  2459  cps.  (C2) 

Differentiating  Af  with  respect  to  time,  we  get 

v2r 

Af  =  — “ —  cos  8  (cos  cos  HR  +  cos  LT  cos  H^.  (C3) 


The  numerical  result  is 


Af  =  0.1381  cpsps. 


First  SAf  is  calculated  using  the  equation 


(C4) 


(C5) 


All  the  quantities  on  the  right-hand  side  are  known,  except  Sf  ,  SL  ,  and  SH  •  For  St  we 
use  a  value  of  lO-10  fT  ,  on  the  basis  of  assumptions  as  to  future  advances  inRelectronicT 
techniques  and  equipment.*  The  errors  in  latitude  and  hour  angle  are  established  on  the 
basis  of  whatever  desired  navigation  accuracy  is  specified.  The  values  used  here  are 


ef  =  io_10f 

T  1 


g  =0.8  naut  mi  x 
lr  R 


S„  =1.3  naut  mi  x 


1852 
R 

1852 

R  cos  1^ 


(C6) 

(C7) 

(C8) 


NRL  Memo  Report  745,  “NRL  Research  and  Development  Program  for  the  Fleet  Ballistic 
Missile:  Summary  of  NRL  Activity  During  August  and  September  1957,”  pp.  23  et  seq.. 
Sept.  30,  1957 


13 


14 


NAVAL  RESEARCH  LABORATORY 


where  R  is  the  earth’s  equatorial  radius  in  meters  and  the  figure  1852  is  the  number  of 
meters  in  one  nautical  mile.  Using  a  transmitted  frequency  fT  of  2200  x  106  cps  and 
r  =  6.378  x  106  meters,  we  get 


=  0.22  cps 

(C9) 

=  2-323  x  10“4  radians 

(CIO) 

=  1.104  x  10”3  radians. 

(Cll) 

The  other  errors  are  given  in  Appendix  A. 

We  now  proceed  to  list  each  term  on  the  right-hand  side  of  Eq.  (C5),  together  with  its 
numerical  value: 

Sf  ^  =  cos  8  ^oos  1^  sin  HR  +  cos  Lj.  sin 


6.048  x  10  14 


_  2 


_  fT^ 


cos  8 


^cos  sin  Hr  +  cos  LT  sin  FL^  8C 


(C12) 


=  1.682  x  10  5 


=  1.486  x  10  3 


(C13) 


(C14) 


(C15) 


fTa>P 


sin  8  ^cos  sin  HR  +  cos  L.J.  sin 


=  1.326  x  10  7 


(Cl  6) 
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'3Af  v  r  fTa,R  /  \  t 

3^  eV  =  L"  cos  8  (sin  **  sin  IIr) \\ 


4.729  x  10-1 


(ir'  s»„)  =  [~^ cos  5  (co* L* cos  “») «»,] 


4.568  x  1(T2 


/  3Af  V  _  ffT"R  /  \  n  "I 

(  -  p  I  -  -  cos  8  (  -  sm  L_  sm  ILr  \  c.T 

vlt  Slt/  c  V  T  V  ltJ 


=  3.029  x  10-7 


( ^  8ht)  =  C°S  8  (C°S  ^  C°S  "0  8ht] 


Adding  these  terms,  we  get 


4.358  x  10”<>. 


(eAf)2  =  0.6338  (cps)2 


and 


®Af  “  0.796  cps. 

Now  S^f  is  calculated  using  the  equation 


(C17) 


(C18) 


(C19) 


(C20) 

(C21) 

(C22) 


it  it  Cj  it  t  it  *)’  t 


(C23) 


Proceeding  as  before, 
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/ a&f  [  fT"2R  f  \  r,~|2 

ecJ  -  j-  — cos  S^cos  Lr  cos  Hr  +  cos  L,.  cos  R..J  £'c 


5.308  x  10-14 


(C25) 


(™  S 

\  ?*>  Co>- 


V  f  2  fTofi  /  \ 

A  -  - - —  cos  $  l  cos  1^  cos  Hr  +  cos  cos  HjJ 


=  1.435  x  10~9 


(Ytf-F  cos  5  ^cos  Lj^  cos  Hr  +  cos  1^  cos 


-|2 


=  4.691  x  10-12 

pC^R 


/  3Af  V  f  fT^R  /  \  l2 

I  —  gg  J  =  -  -  sin  8  f  cos  Lg  cos  HR  +  cos  Lj.  cos  Hjj  gg 


=  4.185  x  10” 1  ** 


a  o  ■ 


8^sin  Lg  cos  ^  gL> 


=  8.118  x  10"11 


(%  =  ["  cos  ?  (*os  ^ sin  Hr)s«r 


=  7.525  x  10“9 


‘ '  [-  ^  ~s  5  (•*"  ^  “■  <)  ^ 


=  3.676  x  10~ls 


CL« 


\2  r  fT6J2R  .  v 

hJ  = - —  co*  s(  cos  lt  sin  “rj  eH, 


(C26) 


(C27) 


(C28) 


(C29) 


(C30) 


(C31) 


=  1.015  x  10-14 


(C32) 


NAVAL  RESEARCH  LABORATORY 


1? 


Adding  these  terms,  we  get 

(^Af)2  =  9.046  x  10-9  (cpsps)  ^ 


and 


~  9.511  x  10  5  cpsps. 


(C33) 

(C34) 


Appendix  D 


FORMULAS  AND  CURVES  FOR  8Af  AND  EA- 

AS  FUNCTIONS  OF  8,  AND  8„ 
lr  hr 

In  order  to  find  out  how  the  requirements  imposed  on  the  accuracy  of  frequency 
measurement  and  frequency  rate  measurement  vary  as  different  prescribed  latitude  and 
longitude  accuracies  are  assumed,  we  take  the  equations  for  SAf  and  gA*f  as  stated  in 
Appendix  B,  hold  everything  at  its  calculated  value  except  the  pertinent  independent  vari¬ 
able  (either  £  or  8H  ),  and  plot  the  four  resulting  curves. 

^R  r 

To  be  more  specific:  if  we  hold  8H^  constant  and  allow  8L^  to  vary,  we  get 

gAf  =  (o.  161  +  8.764  x  106  (Dl) 

and  if  we  hold  8  constant  and  allow  SH  to  vary,  we  get 

R  R 

8Af  =  (o.588  +  3.748  x  104  S^2^172- 
Similarly,  for  gA-f  we  get 

SA'f  3  (s.965  X  10-9  +  1.504  x  10~3  8Lr  ^ 72 

and 

eA‘f  =  (1521  x  10”9  +  6-174  x  10~3  Shr2)1/2  •  (D4) 

The  resulting  curves  are  shown  in  Figs.  Dl  through  D4.  It  is  evident  from  the  curves 
that  a  large  relaxation  in  required  navigation  accuracy  would  give  relatively  small  increases 
in  allowable  doppler  and  doppler  rate  accuracy. 

The  “x”  on  each  curve  indicates  the  original  specified  navigation  error. 


(D2) 

(D3) 
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F  ig,  D1  -  Probable  error  in  doppler  frequency  vs 

prescribed  probable  error  8_  in  latitude 

lr 


Fig.  D2  -  Probable  error  in  doppler  frequency  vs  prescribed 

probable  error  £  in  longitude 

hr 
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Fig,  D3  -  Probable  error  in  doppler  rate  vs 

prescribed  probable  error  0T  in  latitude 


Fig.  D4  -  Probable  error  in  doppler  rate  vs  prescribed 

probable  error  £L  in  longitude 
hr 


Appendix  E 


RESIDUAL  ERROR  IN  H*  AND  L*  ASSUMING  8Af 
AND  gA-f  EQUAL  TO  ZERO 


The  following  notation  is  used  for  the  sake  of  conciseness: 


p  =  tan  Hr 

a  =  a)  c  Af  -  fT  w2!?  cos  8  cos  LT  sin  Hj. 
b  =  c  Af  -  fToj2R  cos  8  cos  LT  cos  HT  . 

Hence 


If  qi  is  any  independent  variable,  then 


(El) 

(E2) 

(E3) 

(E4) 


(E5) 


Assuming  perfect  measurements;  i.e.,  SAf  and  SA*f  equal  to  zero,  there  will  remain 
seven  independent  variables  to  be  considered:  fT  ,  c ,  0>9  r  ,  g,  LJ9  and  Hj..  Therefore  n  =  7 
in  Eq.  (E5). 

To  find  the  required  partial  derivatives,  the  relation 

_  3(arctan  p)  _  1 

1  +  p2  ^qi  {  } 

is  used.  It  is  necessary  first  to  evaluate  p  =  tan  HR ,  so  that  the  common  factor  1/(1  +  p2) 
may  be  known.  If  the  data  in  Appendix  A  and  the  previously  computed  values  of  Af  and  Af 
are  used,  the  result  is 

p  =  -  5.55906.  (E7)* 

Then 

rr; 2  =  3.134x10-2.  (E8) 


The  quantity  p  =  tan  HR  can  also  be  obtained  by  using  sin  HR  and  cos  HR  from  Appendix  A. 
The  values  are  only  slightly  different. 
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Now 


from  Eq.  (E4).  The  numerical  results  for  the  seven  required  terms  are 


Bp  0 

—  =  3.046  x  10~8 

Bfx 

(E10) 

g/Y  =  4.492  x  lO-1? 

\^fT  y 

(Ell) 

Bp 

-  =  -  2.236  X  10  7 

Be 

(E12) 

( —  8  ^  =  1.249  x  10“8 
\  Be  CJ 

(E13) 

Bp 

—  =  1.664  x  106 

B  co 

(E14) 

(f;  8-)2  =  2-769  xl(r4 

(E15) 

Bp 

—  =  1.051  x  10  5 

BR 

(E16) 

(I  "  um  * 10-6 

(E17) 

Bp  i 

—  =  -  2.481  x  101 

BS 

(E18) 

(iF8s)2  =  9-851  x  10-11 

(E19) 

Bp 

-i-  =  -  5.339  x  101 

blt 

(E20) 

=  7.125  * 

(E21) 
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3p  . 

-  =  -  2.887  x  101 

BHj. 


=  8.332  x  10" 10. 


Summing  the  squared  terms  gives 


—  cq-\  -  2.780  X  10'4- 


(E22) 

(E23) 

(E24) 


Extracting  the  square  root  of  Eq.  (E24)  and  multiplying  by  the  value  of  1/(1  +  p2)  from 
Eq.  (E8)  yields 


^  -  5.227  x  10  4  radians, 

where  the  zero  subscript  refers  to  the  fact  that  this  is  the  residual  error, 
is  well  within  the  prescribed  maximum  probable  error 

8„  =  1.104  x  10~3  radians. 

hr 

In  nautical  miles 

^  =  0.616  naut  mi 


(E25) 


This  value 


(Cll) 

(E26) 


compared  with  (see  Appendix  A) 


Consider  now  the  equation 


8„  -  1.3  naut  mi. 

hr 


cos  Lr  =  ( 1  +  tan2 


(l  +  tan2  hY/2  fcos  LT  cos  ^ - 

\  *v  fjco2  R  cos  S  / 


It  is  convenient  to  introduce  the  notation 

r  =  cos  L„ 


s  -  cos  LT  cos  Hj. 


c  Af 


Then  Eq.  (25)  can  be  written 


fT  a)2  R  cos  8 


»(1  +  P2) 


1/2 


(E27) 

(25) 

(E28) 

(E29) 

(E30) 


If  q.  is  any  of  the  independent  variables,  then 

BLR  =  9(arccos  r)  _  _  f2  -1/2  Br 

9q.  9q.  3q. 


(E31) 
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It  is  necessary  to  calculate  r  =  s(l  +  p2)1/2  and  from  that  to  get  (1  -  r2)  1/2  . 

The  results  are 

s  =  6.05421  X  10-2  (E32)* 

(1  +  p2)1/2  =  5.64829.  (E33) 

Hence 

r  =  3-41959  x  10_1  (E34) 

(1  -  r 2)-l/2  =  1.06417.  (E35) 

From  Eq.  (E30) 


Br 

Bqi 


=  a  +  p2)1/2  -g 


+ 


sp(l  +  P2) 


-1/2  _l£L 

da.  ' 


(E36) 


The  only  quantity  which  one  has  to  compute  is  Bs/Bq. ,  since  all  the  other  quantities  are 
known;  Bp/Bq.  and  p  are  obtained  from  the  above  hr  calculations.  Putting  in  the  numerical 
values,  we  get 

=  5.64829  -5.95858  x  10”2  •  (E37) 

3q.  3q.  3q.  ' 

Evaluation  of  the  seven  terms  yields 


=  -2.961  x  10_1° 

(E38) 

=  4.244  x  10-21 

(E39) 

=  2.173  x  10-9 

(E40) 

=  1.180  x  10~12 

(E41) 

=  -7.499  x  103 

(E42) 

=  5.623  x  10"9 

(E43) 

^Strictly  speaking,  (1  +  p2)1/2  =  sec  HR,  which  is  negative.  This  introduces  a  negative 
sign  into  both  terms  of  Eq.  (E36).  Since  each  Br/Bqi  is  to  be  squared,  the  positive 
(1  +  p2)1/2  was  used. 
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Br 

—  =  -1.022  x  10  7 
BR 


&')2- 


1.044  x  10~10 


—  =  2.412  x  10  1 

38 


309  x  10"15 


Br 

—  =  2.466  x  10  1 
BL_ 


- 8t  1  =  1.520  x  10  14 


Br 

BRj. 


f*.  '2 


=  -7.183  x  10'1 


V3HT 

Summation  of  the  squared  terms  gives 


gj,  \  =  5.159  x  10-13  . 


7.236  x  10  9 . 


Extraction  of  the  square  root  and  multiplication  by  (i  -  r2)-1/2  =  1.06417  gives 

(£.  \  =  8.506  x  10-5  radians 

V  */° 

in  absolute  value.  This  may  be  compared  with  the  prescribed  value 


(E44) 

(E45) 

(E46) 

(E47) 

(E48) 

(E49) 

(E50) 

(E51) 

(E52) 

(E53) 


In  nautical  miles 


=  2.323  x  10  4  radians. 


=  0.293  naut  mi 


compared  with  (see  Appendix  A) 


p  =0.8  naut  mi. 
lr 
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